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Abstract 

Background: Organisms may respond to environmental change by means of genetic adaptation, phenotypic 
plasticity or both, which may result in genotype-environment interactions (G x E) if genotypes differ in their 
phenotypic response. We here specifically target the latter source of variation (i.e. G x E) by comparing plastic 
responses among lines of the tropical butterfly Bicyclus anynana that had been selected for increased cold 
tolerance and according controls. Our main aim here was to test the hypothesis that directional selection on cold 
tolerance will interfere with plastic capacities. 

Results: Plastic responses to temperature and feeding treatments were strong, with e.g. higher compared to lower 
temperatures reducing cold tolerance, longevity, pupal mass, and development time. We report a number of 
statistically significant genotype-environment interactions (i.e. interactions between selection regime and 
environmental variables), but most of these were not consistent across treatment groups. We found some evidence 
though for larger plastic responses to different rearing temperatures in the selection compared to the control lines, 
while plastic responses to different adult temperatures and feeding treatments were overall very similar across 
selection regimes. 

Conclusion: Our results indicate that plastic capacities are not always constrained by directional selection (on cold 
tolerance) and therefore genetic changes in trait means, but may operate independently. 

Keywords: Artificial selection, Bicyclus anynana, Constraint, Genetic adaptation. Genotype by environment interaction, 
Phenotypic plasticity. Temperature stress resistance 



Background 

Temperature is considered one of the most important 
selective agents, and consequently research on tem- 
perature stress resistance has attracted much interest 
over recent decades [1-5]. Temperature-stress resistance 
refers to an organism's ability to cope with stressfully 
high or low temperatures, and is considered a key factor 
for explaining the distribution and abundance of species 
[6,7]. Enhanced resistance to temperature stress can 
be reached by means of phenotypic plasticity, i.e. non- 
genetic physiological changes as a direct response to 
environmental variation, or genetic adaptation [5,8]. 
Genetic variation in temperature stress resistance has 
been commonly reported in natural systems, with e.g. 
tropical species showing a higher heat but a lower cold 
tolerance than temperate-zone species and vice versa 
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[1,3,4,6,7,9-13]. Likewise, high-altitude populations typic- 
ally show a lower heat but a higher cold tolerance than 
low- altitude populations [14]. Such geographic variation 
in fitness-related traits provides strong evidence that 
these patterns have been shaped by natural selection 
[14,15]. Several species are also known to respond readily 
to artificial selection on thermal tolerance traits, provid- 
ing direct experimental proof for genetic adaptation in 
temperature stress resistance (e.g. [2,16,17]). 

In addition to genetic adaptation, phenotypic plasticity 
provides a further efficient mechanism to cope with 
temperature variation [18]. We here define phenotypic 
plasticity as the set of different phenotypes that may be 
produced by a single genotype in direct response to differ- 
ent environmental conditions [8]. Plastic responses can be 
induced during development (developmental plasticity) or 
in the adult stage (adult acclimation; [5,19,20]), and bene- 
ficial effects of a brief exposure to less extreme tempera- 
tures are referred to as rapid hardening [1,21,22]. Studies 
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on plastic responses in temperature stress resistance 
revealed that cooler rearing or adult temperatures increase 
cold but decrease heat resistance and vice versa, compris- 
ing a near universal pattern of putatively adaptive pheno- 
typic plasticity [1,6,23]. For instance, groups of a tropical 
butterfly reared in a common environment, but being 
exposed for two days to either 20°C or 27°C in the adult 
stage, differed subsequently by a factor 2-3 in heat toler- 
ance [23]. Thus, plastic changes in thermal tolerance may 
be induced within short periods of time, and are highly ef- 
fective means to deal with temperature variation. 

Both genetic variation and phenotypic plasticity may 
finally interact with one another, resulting in genotype- 
environment interactions (G x E). We here use G x E to 
refer to different degrees by which individual genotypes 
are able to respond to environmental variation [24]. 
Consequently, a significant G x E demonstrates that 
some genotypes are more plastic than others [14]. Quan- 
tifying such variation is consequently used to explore 
genetic variation in phenotypic plasticity [25-28]. G x E 
has been reported for e.g. insect growth and develop- 
ment, fecundity, and immune function [29,30], and is 
considered to contribute substantially to the mainten- 
ance of the genetic variation of traits [31,32]. Although 
interactions between genotypes and environmental con- 
ditions may thus represent a significant source of vari- 
ation, our respective knowledge is still fairly limited. 

Against this background we here investigate G x E in 
control and selection lines of the tropical nymphalid 
butterfly Bicyclus anynana (Butler, 1879). This species is 
known to respond readily to temperature manipulations 
by phenotypic plasticity in temperature stress resistance 
as well as in an array of other traits [20,23]. Different 
genotypes of B. anynana had been previously produced 
by applying artificial selection to chill-coma recovery 
time, yielding highly divergent lines [17]. Chill-coma re- 
covery time is the time an individual needs to regain 
mobility after cold exposure, and is considered a reliable 
proxy of climatic cold adaptation [6,23,32,33]. This trait 
is furthermore known to respond readily to selection as 
well as to ambient temperatures [2,16,23,34]. We thus 
exploit the genetic variation generated through artificial 
selection to investigate whether genotypes substantially 
varying in cold tolerance respond differentially to envir- 
onmental manipulations. This explicit focus on G x E 
contrasts with earlier studies from our group investigat- 
ing responses and correlated responses to selection [34], 
effects of selection in the adult stage throughout the life 
cycle [35], environmental effects on cold tolerance [19], 
or effects of inbreeding [17]. Specifically we here test the 
hypotheses that lines selected for increased cold toler- 
ance show reduced plastic responses in cold tolerance 
(and possibly further traits) compared to control lines. 
This expectation rests on the assumption that the cold- 



tolerant lines should benefit less from strong plastic 
responses in case of critically low temperatures com- 
pared to control lines. Further, trait values may already 
be close to their physiologically determined maximum, 
constraining further plastic increases. 

In order to induce plastic responses we exposed selec- 
tion and control lines to different rearing and acclima- 
tion temperatures as well as to different feeding regimes, 
and measured the respective effects on cold tolerance 
and additionally on development and longevity. Some of 
the concomitant effects reported here are hardly novel, 
though including environmental manipulations is evi- 
dently necessary to investigate G x E. Furthermore this 
enables us to test a number of additional hypotheses, 
namely that higher temperatures are predicted to speed 
up development and to reduce pupal mass, longevity, 
and cold tolerance [5,20,36,37]. A period of starvation is, 
based on earlier results, predicted to reduce longevity 
but to leave cold tolerance unaffected ([23,36,37]; see 
Table 1 for specific hypotheses). 

As the initial design of the selection experiment 
included inbreeding as a factor, we furthermore consider 
effects of inbreeding in our study. However, based on the 
fact that the lines were allowed to mate randomly after the 
different inbreeding levels had been established, no effects 
of inbreeding are expected anymore (see further below). 
Selection regime is, based on previous results, also 
expected to yield no effects on traits other than cold toler- 
ance [34]. We finally include sex as a factor in our ana- 
lyses, predicting that females show a slower larval but a 
faster pupal development, a higher pupal mass, and a 
higher longevity compared to males ([23], Table 1). 

Methods 

Study organism 

Bicyclus anynana (Butler 1897; Nymphalidae, Satyrinae) 
is a tropical fruit-feeding butterfly, ranging from southern 
Africa to Ethiopia [38]. The species inhabits regions with 
alternating dry and wet seasons, and shows accordingly 
two seasonal morphs. During the colder dry season (ca. 
18°C, May until November) the species has rather uniform 
wing patterns and small eyespots, while it exhibits large 
eyespots and bright bands on both wings in the warmer 
wet season (ca. 23°C, December until April; [39]). During 
the dry season reproduction ceases and butterflies do not 
mate before the onset of the next wet season [29] . Females 
are relatively monandrous, though multiple mating occurs 
in the field and in the laboratory [40]. A stock population 
of B. anynana was founded at Greifswald University in 
2007 from several hundred individuals derived from a 
well-established stock population at Leiden University, 
The Netherlands. The latter was founded in 1988 from 
over 80 gravid females collected at a single locality in 
Nkhata Bay, Malawi. To maintain high levels of 
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Table 1 Hypotheses for the main effects of several factors on larval time, pupal time, larval growth rate, pupal mass, 
chill-coma recovery time, and longevity 



Traits 






Factors 








Selection regime 


Inbreeding 


Rearing temperature 


Sex 


Adult temperature 


Feeding regime 


Larval time 


No 


No 


20°C > 27°C 


Females > Males 






Pupal time 


No 


No 


20°C > 27°C 


Males > Females 






Growth rate 


No 


No 


27°C > 20°C 


Males > Females 






Pupal mass 


No 


No 


20°C > 27°C 


Females > Males 






Chill-coma recovery 
time 


Controls > Selected Lines 


No 


27°C > 20°C 


No 


27°C > 20°C 


No 


Longevity 


No 


No 


20°C > 27°C 


Females > Males 


20°C > 27°C 


Control > Starvation 



The hypotheses stated rest on prior results mainly obtained from B. anynana. Selection for increased cold tolerance is expected to have no effect on any trait 
other than chill-coma recovery time, as no correlated responses to selection were found in a previous study [34]. No inbreeding effects are expected here due to 
random mating for 14 generations after full-sib matings and thus prior to scoring trait values. Hypotheses regarding temperature, sex and feeding effects are all 
straight-forward and rest on previous results [e.g. 5,17,20,23,34,36,37]. '-': Not tested. 



heterozygosity several hundred adults are reared in each 
generation [41]. For this experiment animals from the 
Greifswald stock population were used. 

Experimental design 

We here used 12 selection lines that had been previously 
established at the Department of Animal Ecology, Greifs- 
wald University [34,35]. Selection lines for increased cold 
stress resistance (shorter chill-coma recovery times) and 
according unselected controls had been derived from 
three different levels of inbreeding, using a full-sib breed- 
ing design: outbred controls (C) having resulted from mat- 
ings between unrelated butterflies, inbreeding 1 (II) 
having resulted from matings between full sibs, and 
inbreeding 2 (12) having resulted from matings between 
full sibs in two consecutive generations. We used ca. 120 
full-sib families each for inbreeding levels C, II, and 12 
[34]. At the start of the selection experiment though 
butterflies were pooled across families within inbreeding 
levels. Thus, during the course of the selection experiment 
(10 generations) mating was random within lines. Per 
inbreeding level, four lines were set up, two for increased 
cold stress resistance and two unselected controls. This 
design resulted in a total of 6 selection and 6 control lines 
(for details see [17,34]). Lines had been kept without selec- 
tion under standard rearing conditions for 4 generations 
prior to this experiment. Several hundred butterflies were 
reared per line in each generation. For the current experi- 
ment eggs were collected from all 12 lines. Eggs were 
thereafter randomly divided among a low (20°C) and a 
high (27°C) rearing temperature (70 ± 5% relative humid- 
ity and photoperiod of L12:D12 throughout). The two 
temperatures chosen are similar to the ones this species 
experiences during the wet and dry season in the field, re- 
spectively [29]. Larvae were reared in sleeve-like gauze 
cages, using 10 replicate cages per line and rearing 
temperature (resulting in a total of 240 cages), and a 
standard density of 20 (20°C) and 40 (27°C) larvae per 



cage, respectively. The higher density per cage at 27°C was 
due to the need for more individuals at this temperature 
(see below). Care was taken though that developing larvae 
were never exposed to any food shortage, being fed on 
young maize plants ad libitum throughout. Resulting 
pupae were collected daily, weighed one day after pupa- 
tion, and were then individually transferred to small plas- 
tic cups (volume 125 ml). For all animals we scored larval 
time (from egg-laying until pupation, thus including egg 
development), pupal time (from pupation until adult eclo- 
sion), larval growth rate (In pupal mass/larval time), and 
pupal mass (measured on day 2 after pupation). 

Following adult eclosion, all butterflies were marked 
individually and afterwards once again randomly divided 
among 20°C and 27°C, resulting in four rearing by adult 
temperature groups per selection line (20-20°C, 20-27°C, 
27-20°C, 27-27°C). While the butterflies reared at 27°C 
were a last time divided among two feeding treatments, 
being fed with banana (control) or water only (starva- 
tion), all animals reared at 20°C were fed with banana ad 
libitum (to keep the size of the experiment manageable). 
On day two after eclosion, we tested 24-52 individuals 
per sex, treatment group and line for chill-coma recov- 
ery time, resulting in a total number of 4772 butterflies. 
Therefore, all butterflies were transferred individually to 
plastic cups (125 ml), arranged on a tray in a rando- 
mized block design, and then exposed to 1°C for 19 h to 
induce a chill coma. Recovery time (the time until the 
butterflies were able to stand up) was scored in a climate 
cell at 20°C. Observations were terminated after 60 minutes. 
Butterflies that had not recovered by then were given 
the maximal recovery time of 60 minutes. This applied 
to < 5% of the butterflies tested. 

Thus, we used the same technique for measuring chill- 
coma recovery time as during selection [35]. This trait is 
considered a reliable proxy of climatic cold adaptation, 
is largely independent of the method used to induce a 
chill coma, and has been used successfully in B. anynana 
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before [24]. To afterwards score longevity the butterflies 
reared at 20°C were returned to their respective adult 
temperatures and fed ad libitum, while the butterflies 
reared at 27°C were all transferred to 27°C (due to space 
limitations) and fed as outlined above until death. 

Statistical analyses 

Data on larval time, pupal time, pupal mass, growth rate, 
chiU coma recovery time, and longevity were analyzed 
with mixed model analyses of variance (ANOVAs), with 
selection regime (selection vs. control), level of inbreeding 
(C, II, 12), rearing temperature (20 vs. 27°C), adult 
temperature (20 vs. 27°C), adult feeding regime (control 
vs. starvation, for animals reared at 27°C only), and/or sex 
as fixed factors. Replicate line was nested within selection 
regime and inbreeding level, and replicate cage was nested 
within replicate line, inbreeding level and selection regime. 
Both latter factors were included as random effects. Mini- 
mum adequate models were constructed by sequentially 
removing non-significant interaction terms. Due to differ- 
ences in the experimental set-up across rearing tempera- 
tures (see above), three different analyses were run for 
chiU-coma recovery time and longevity. For chill-coma re- 
covery time we analyzed all individuals together (note that 
effects of adult feeding regime were non-significant) and 
separated by rearing temperature, as different adult feed- 
ing regimes were only applied to animals reared at 27°C. 
For longevity, effects of rearing temperature were analyzed 
in animals fed ad libitum and kept as adults at 27°C, 
effects of adult temperature were analyzed in animals 
reared at 20°C (as only in these different adult tempera- 
tures were used throughout), and effects of adult feeding 
regime were analyzed in animals reared at 27°C (as only in 
these different adult feeding regimes were applied). Sig- 
nificant differences between inbreeding levels were located 
using Tukeys HSD posthoc test. When presenting the 
statistical results we start with the main effects and con- 
tinue with interactions, which can be problematic in case 
of strong interactive effects. Therefore, main effects are 
interpreted cautiously in case of significant interactions. 
All statistical tests were performed by using JMP (4.0.0) or 
Statistica (6.1). Throughout, all means are given ± 1 SE. 

Results 

Larval time 

Larval time was significantly affected by rearing temperature, 
sex, replicate line and replicate cage, but not by selection 
regime and inbreeding level (Table 2A). Larval time was 
much longer at the lower compared to the higher rearing 
temperature (50.6 ± 0.06 d > 27.4 ± 0.05 d; Figure lA), 
and was longer in females than in males (40.2 ± 0.06 d > 
37.8 ± 0.05 d). The absolute difference in larval time be- 
tween males and females was smaller at 27°C (26.4 ± 0.06 
d vs. 28.3 ± 0.07 d) than at 20°C (49.2 ± 0.09 d vs. 52.0 ± 



0.09 d), as indicated by the significant rearing temperature 
by sex interaction. The sexual difference was particularly 
small in the unselected controls reared at 27°C (significant 
selection regime * rearing temperature * sex interaction; 
Figure IB). The significant interaction between selection 
regime and rearing temperature reveals a significantly 
longer larval time in control than in selection lines at 20°C 
(51.0 ± 0.09 d > 50.2 ± 0.09 d), which was not the case at 
27°C (27.5 ± 0.06 d = 27.3 ± 0.07 d; Tukey HSD). At the 
higher rearing temperature, larval time was longer in the 
control compared to the inbred groups (C: 28.0 ± 0.08 d > 
II: 27.1 ± 0.08 d = 12: 27.1 ± 0.08 d), while there were no 
significant differences at 20°C (C: 50.5 ± 0.1 d = II: 50.3 ± 
0.1 d = 12: 50.9 ± 0.1 d; Tukey HSD; significant inbreeding 
level * rearing temperature interaction). Both above pat- 
terns though were not consistent across all groups (signifi- 
cant selection regime * inbreeding level * rearing 
temperature interaction; Figure lA). 

Pupal time 

Pupal time was significantly influenced by rearing 
temperature, sex, replicate cage and replicate line, but not 
by selection regime and inbreeding level (Table 2B). Pupal 
time was longer at the lower compared to the higher rear- 
ing temperature (11.9 ± 0.02 d > 7.0 ± 0.01 d), and was 
longer in males than in females (9.7 ± 0.02 d > 9.2 ± 0.02 
d). A significant rearing temperature by sex interaction 
indicates that sexual differences were slightly more pro- 
nounced at 27''C (7.4 ± 0.02 d > 6.7 ± 0.02 d) than at 20°C 
(12.1 ± 0.03 d > 11.6 ± 0.03 d; Figure IC). At the higher 
rearing temperature, pupal time increased with increasing 
inbreeding level (C: 6.8 ± 0.02 d < II: 7.1 ± 0.02 d < 12: 7.2 
± 0.02 d), while at the lower rearing temperature inbreed- 
ing level had no significant effect (C: 11.9 ± 0.03 d = II: 
11.9 ± 0.03 d = 12: 11.8 ± 0.03 d; Tukey HSD; significant 
inbreeding level * rearing temperature interaction). Fur- 
thermore, inbreeding effects on pupal time were more 
pronounced in males (C: 9.6 ± 0.03 d < II: 9.8 ± 0.03 d = 
12: 9.8 ± 0.03 d) than in females (C: 9.1 ± 0.03 d = II: 9.2 
± 0.03 d = 12: 9.2 ± 0.03 d; Tukey HSD; significant 
inbreeding level * sex interaction). Consequently, inbreed- 
ing effects were largely restricted to males reared at 27°C 
(significant inbreeding level * rearing temperature * sex 
interaction; Figure IC). The significant interaction be- 
tween selection regime, inbreeding level and rearing 
temperature indicates that differences between selection 
and control lines were basically limited to the outbred 
control group when reared at 27°C (Figure ID). 

Growth rate 

Growth rate was also significantly affected by rearing 
temperature, sex, replicate line and replicate cage, but not 
by selection regime and inbreeding level (Table 2C). 
Growth rate was higher at the higher compared to the 
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Table 2 Results of nested ANOVAs for larval time (A), 
pupal time (B), larval growth rate (C) and pupal mass (D) 



A) Larval time 


DF 


MS 


F 


P 


Selection Regime 


1,6 


205.6 


1.0 


0.3453 


hbreeding Level 


2,6 


87.2 


04 


0.6614 


Replicate Line [Sel. Reg. 

& Inbr.] 


6,110 


216.0 


94 


<0.0001 


Repl. Cage [Sel., Inbr. 
& Repl.] 


108,4640 


23.5 


3.8 


<0.0001 


Rearing Temperature 


1,4640 


555906.0 


90265.2 


<0.0001 


Sex 


1,4640 


5780.5 


922.0 


<0.0001 


Selection Reg. * Inbreeding 
Level 


2,6 


236.5 


1.2 


0.3635 


Selection Reg. * Rearing Temp. 


1,4640 


1074 


17.1 


<0.0001 


Selection Reg. * Sex 


1,4640 


11.2 


1.8 


0.1819 


Inbreeding Level * Rearing 
Temp. 


24640 


139.2 


22.2 


<0.0001 


nbreeding Level * Sex 


24640 


144 


2.3 


0.1007 


Rearing Temp. * Sex 


1,4640 


219.6 


35.0 


<0.0001 


Sel. Reg. * Inbreed. * Rear. 
Temp. 


2,4640 


49.9 


8.0 


0.0004 


Sel. Reg. * Rearing Temp. 
* Sex 


1,4640 


25.0 


25.0 


0.0457 


Error 


4640 


63 






B) Pupal time 


DF 


MS 


F 


P 


Selection Regime 


1,6 


14 


01 


0.7627 


nbreeding Level 


2,6 


9.2 


0.6 


0.5567 


Replicate Line [Se . Reg. 
& Inbr.] 


6,111 


15.6 


1 1.7 


<0.0001 


Repl. Cage [Sel., Inbr. 
& RepL] 


108,4639 


14 


2.5 


<0.0001 


Rearing Temperature 


1,4639 


24740.3 


45323.1 


<0.0001 


Sex 


1,4639 


3764 


689.6 


<0.0001 


Selection Reg. * Inbreeding 
Level 


2,6 


1.7 


01 


0.8888 


Selection Reg. * Rearing 
Temp. 


1,4639 


0.2 


04 


0.5196 


Selection Reg. * Sex 


1,4639 


0.8 


8.8 


0.3708 


hbreeding Level * Rearing 
Temp. 


24639 


12.8 


234 


<0.0001 


hbreeding Level * Sex 


24639 


3.0 


5.5 


0.0041 


Rearing Temp. * Sex 


1,4639 


13.6 


24.9 


<0.0001 


Sel. Reg. * Inbreed. * Rear. 
Temp. 


2,3941 


3.5 


6.3 


0.0018 


hbreed. * Rear. Temp. 
* Sex 


24639 


2.2 


4.1 


0.0173 


Error 


4639 


0.5 






C) Growth rate 


DF 


MS 


F 


P 


Selection Regime 


1,6 


5.7 


01 


0.7710 


hbreeding Level 


2,6 


40.9 


0.7 


0.5514 



Table 2 Results of nested ANOVAs for larval time (A), 
pupal time (B), larval growth rate (C) and pupal mass (D) 

(Continued) 



Replicate Line [Sel. Reg. 
& Inbr.] 


6,109 


68.1 


8.2 


<0.0001 


Repl. Cage [Sel., Inbr. 
& Repl.] 


108,4641 


8.6 


4.3 


<0.0001 


Rearing Temperature 


1,4641 


79417.9 


40343.3 


<0.0001 


Sex 


1,4641 


154.2 


78.3 


<0.0001 


Selection Reg. * Inbreeding 
Level 


2,6 


43.2 


0.7 


0.5353 


Selection Reg. * Rearing 
Temp. 


1,4641 


2.1 


1.1 


0.3012 


Selection Reg. * Sex 


1,4641 


0.4 


0.2 


0.6716 


Inbreeding Level * Rearing 
Temp. 


2,4641 


33.0 


16.8 


<0.0001 


Inbreeding Leve * Sex 


2,4641 


0.6 


0.3 


0.7479 


Rearing Temp. * Sex 


1,4641 


31.3 


15.9 


<0.0001 


Sel. Reg. * Inbreed. * Rear. 
Temp. 


24641 


8.8 


4.5 


0.0115 


Error 


4641 


2.0 






D) Pupal mass 


DF 


MS 


F 


P 


Selection Regime 


1,6 


2366.2 


0.6 


0.4880 


Inbreeding Level 


2,6 


38734.6 


8.9 


0.0156 


Replicate Line [Sel. Reg. 
& Inbr.] 


6,110 


4746.8 


3.2 


0.0064 


Repl. Cage [Sel., Inbr. 
& Repl.] 


108,4641 


15274 


3.8 


<0.0001 


Rearing Temperature 


1,4641 


1 56622.0 


389.2 


<0.0001 


Sex 


1,4641 


1661991.0 


4143.6 


<0.0001 


Selection Reg. * Inbreeding 
Level 


2,6 


2314 


0.1 


0.9487 


Selection Reg. * Rearing 
Temp. 


1,4641 


949.7 


24 


01239 


Selection Reg. * Sex 


1,4641 


2514 


0.6 


0.4286 


Inbreeding Level * Rearing 
Temp. 


24641 


1285.1 


3.2 


0.0407 


Inbreeding Level * Sex 


24641 


1253.2 


3.1 


0.0441 


Rearing Temp. * Sex 


1,4641 


2954 


0.7 


0.3908 


Sel. Reg. * Inbreed. * Rear. 
Temp. 


24641 


2182.8 


54 


0.0044 


Error 


4641 


400.5 







Results for the effects of selection regime (selection vs. control), inbreeding 
level (C, II, 12}, replicate line, replicate cage, rearing temperature (20 vs. 27°C), 
and sex on larval time (A), pupal time (B), larval growth rate (C), and pupal 
mass (D) in Bicyclus anynana. Replicate line was nested within selection regime 
and inbreeding level, and replicate cage was nested within replicate line, 
inbreeding level and selection regime. Both latter factors were included as 
random effects, while all others were considered fixed effects. Minimum 
adequate models were constructed by sequentially removing non-significant 
interaction terms. Significant P-values are given in bold. 
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Figure 1 Variation in life-history traits (means + 1 SE) in relation to various factors in Bicyclus anynana. Larval time in relation to rearing 
temperature, inbreeding level and selection regime (A), and in relation to rearing temperature, selection regime and sex (B); pupal time in 
relation to rearing temperature, inbreeding level and sex (C), and in relation to rearing temperature, inbreeding level and selection regime (D); 
larval growth rate in relation to rearing temperature, inbreeding level and selection regime (E); and pupal mass in relation to rearing temperature, 
inbreeding level and selection regime (F). UC: unselected control; SL: selection line. Please note that in some cases standard errors are so smal 
that error bars are not visible. 



lower rearing temperature (19.1 ± 0.03%/d > 10.4 ± 
0.03%/d), and was higher in males than in females (14.9 ± 
0.03%/d > 14.6 ± 0.03%/d). The significant interaction be- 
tween rearing temperature and sex indicates a larger sex- 
ual difference in growth rates at 27°C (males: 19.4 ± 
0.04%/d > females: 18.8 ± 0.04%/d) than at 20°C (males: 
10.5 ± 0.05%/d > females: 10.3 ± 0.05%/d). Effects of 
inbreeding were found at the higher rearing temperature 
only (C: 18.8 ± 0.04%/d < II: 19.3 ± 0.04%/d = 12: 19.2 
± 0.05%/d), but not at the lower one (C: 10.4 ± 0.06%/d 
= II: 10.5 ± 0.06%/d = 12: 10.3 ± 0.06%/d; Tukey HSD, 
significant inbreeding level * rearing temperature inter- 
action). Moreover, clear evidence for effects of inbreed- 
ing on growth rates was restricted to the unselected 
control lines reared at 27°C (significant selection regime 
* inbreeding level * rearing temperature interaction; 
Figure IE). 

Pupal mass 

Pupal mass differed significantly between inbreeding 
levels, rearing temperatures, sexes, replicate lines and 



replicate cages, but not in relation to selection regime 
(Table 2D). Pupal mass was highest in the inbreeding 1 
group, followed by the outbred control and finally the 
inbreeding 2 group (II: 192.7 ± 1.72 mg > C: 188.9 ± 
1.74 mg > 12: 182.1 ± 1.85 mg; Tukey HSD). The exact 
patterns though differed among the higher (II: 186.0 ± 
0.63 mg = C: 183.9 ± 0.64 mg > 12: 175.4 ± 0.65 mg) and 
the lower rearing temperature (U: 199.3 ± 0.83 mg > C: 
193.9 ± 0.85 mg > 12: 188.7 ± 0.91 mg; Tukey HSD; signifi- 
cant inbreeding level * rearing temperature interaction), 
and also among the sexes (males: II: 173.1 ± 0.73 mg > C: 
168.0 ± 0.73 mg > 12: 163.0 ± 0.77 mg; females: II: 212.2 ± 
0.74 mg = C: 209.7 ± 0.74 mg > 12: 201.2 ± 0.78 mg; sig- 
nificant inbreeding level * sex interaction). Effects of selec- 
tion regime were absent except for the outbred controls 
reared at 27°C, where pupal mass was higher in the unse- 
lected controls (significant selection regime * inbreeding 
level * rearing temperature interaction; Figure IF). Pupal 
mass was overall higher at 20°C than at 27°C (194.0 ± 
0.50 mg > 181.8 ± 0.37 mg), and was higher in females 
than in males (207.7 ± 0.44 mg > 168.0 ± 0.43 mg). 
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Table 3 Results of nested ANOVAs for chill-coma recovery 
time 



A) Chill-coma recovery time 


DF 


MS 


F 


P 


Selection Regime 


1,6 


24368.4 


22.0 


0.0032 


hbreeding Level 


2,6 


1 1 834 


1.1 


0.4005 


Replicate Line [Sel. & Inbr.] 


6,4738 


1 204.3 


7.7 


<0.0001 


Rearing Temperature 


1,4738 


127708.0 


815.0 


<0.0001 


Adult Temperature 


1,4738 


1 1 5643.0 


738.0 


<0.0001 


Sex 


1,4738 


<0.1 


<0.1 


0.9926 


Selection Reg. * Inbreeding Level 


2,6 


192.6 


0.2 


0.8448 


Selection Reg.* Rearing Temp. 


1,4738 


480.4 


3.1 


0.0800 


Selection Reg. * Adult Temp. 


1,4738 


1327.2 


8.5 


0.0036 


Selection Reg. * Sex 


1,4738 


21.2 


0.1 


0.7131 


nbreeding Level * Rearing Temp. 


2,4738 


797.1 


5.1 


0.0062 


hbreeding Level * Adult Temp. 


2,4738 


6924 


44 


0.0121 


hbreeding Level * Sex 


2,4738 


254.0 


1.6 


0.1979 


Rearing Temp. * Adult Temp. 


1,4738 


1.2 


<0.1 


0.9298 


Rearing Temp. * Sex 


1,4738 


1.8 


<0.1 


0.9140 


Adult Temp. * Sex 


1,4738 


1 384.2 


8.8 


0.0030 


Sel. Reg. * Inbreed. * Rear. Temp. 


2,4738 


910.95 


5.8 


0.0030 


Sel. Reg. * Inbreed. * Adult Temp. 


2,4738 


527.1 


34 


0.0347 


hbreed. * Rear. Temp. * Adult Temp. 


2,4738 


504.9 


3.2 


0.0400 


Rear. Temp. * Adult Temp. * Sex 


1,4738 


1121.7 


7.2 


0.0075 


Error 


4738 


156.7 






B) Animals reared at 20°C 


DF 


MS 


F 


P 


Selection Regime 


1,6 


6782.9 


18.8 


0.0049 


hbreeding Level 


2,6 


213.8 


0.6 


0.5829 


Replicate Line [Sel. & Inbr.] 


6,1668 


362.6 


4.1 


0.0004 


Adult Temperature 


1,1668 


44483.1 


501.1 


<0.0001 


Sex 


1,1668 


1.3 


<0.1 


0.9046 


Selection Reg. * Adult Temp. 


1,1668 


1458.2 


164 


<0.0001 


Error 


1668 


88.8 






C) Animals reared at 27°C 


DF 


MS 


F 


p 


Selection Regime 


1,6 


23033.8 


15.3 


0.0079 


Inbreeding Level 


2,6 


2340.1 


1.5 


0.2859 


Replicate Line [Sel. & Inbreed.] 


6,3059 


1 506.7 


7.9 


<0.0001 


Adult Temperature 


1,3059 


81814.6 


427.7 


<0.0001 


Adult Feeding Regime 


1,3059 


5.6 


<0.1 


0.8644 


Sex 


1,3059 


04 


<0.1 


0.9620 


Selection Reg. * Inbreeding Level 


2,6 


959.2 


0.6 


0.5611 


Selection Reg. * Adult Temp. 


1,3059 


321.2 


1.7 


0.1952 


Selection Reg. * Adult Feed. 


1,3059 


66.0 


0.3 


0.5570 


Selection Reg. * Sex 


1,3059 


38.0 


0.2 


0.6557 


hbreeding Level * Adult Temp. 


2,305 


1 643.0 


8.6 


0.0002 


hbreeding Level * Adult Feed. 


2,3059 


211.2 


1.1 


0.3317 


hbreeding Level * Sex 


2,3059 


145.2 


0.8 


0.4682 



Table 3 Results of nested ANOVAs for chill-coma recovery 
time (Continued) 



Adult Temp. * Adult Feed. 


1,3059 


231.8 


1.2 


0.2711 


Adult Temp. * Sex 


1,3059 


3459.7 


18.1 


<0.0001 


Adult Feed. * Sex 


1,3059 


174.4 


0.9 


0.3397 


Sel. Reg. * Inbreed. * Adult Temp. 


2,3059 


961.3 


5.0 


0.0066 


Inbreed. * Adult Feed. * Sex 


2,3059 


1 2434 


6.5 


0.0015 


Adult Temp. * Adult Feed. * Sex 


1,3059 


1362.7 


7.1 


0.0076 


Error 


3059 


191.3 







Results for the effects of selection regime (selection vs. control), inbreeding 
level (C, II, 12}, replicate line, rearing temperature {20 vs. 27°C), adult 
temperature (20 vs. 27''C), and sex on chill-coma recovery time in Bicyclus 
anynana (A). In Tables B (animals reared at 20°C) and C (animals reared at 
27°C) data were analyzed separately for each rearing temperature, as different 
feeding treatments were employed in animals reared at 27°C only. 
Consequently, the factor rearing temperature was dropped from both 
analyses, and the factor adult feeding regime (control vs. starvation) was 
added to Table C. For details see Material and Methods. Replicate line (random 
factor) was nested within selection regime and inbreeding level throughout. 
Minimum adequate models were constructed by sequentially removing 
non-significant interaction terms. Significant P-values are given in bold. 

Chill-coma recovery time 

Chill-coma recovery time was significantly affected by 
selection regime, rearing temperature, adult tem- 
perature, and replicate line, but not by inbreeding level 
and sex (Table 3A). Animals from the selection 
lines recovered faster than those from the control lines 
(18.2 ± 0.73 min > 23.0 ± 0.70 min). Furthermore, ani- 
mals reared at 20°C recovered faster than those reared at 
27''C (15.2 ± 0.31 min < 26.0 ± 0.23 min), as was the 
case for animals acclimated to 20°C compared to 27°C 
(15.4 ± 0.27 min < 25.8 ± 0.27 min). The significant 
interaction between selection regime and adult 
temperature reveals that the difference between selection 
lines and unselected controls was larger at 27°C (22.9 ± 
0.38 min versus 28.7 ± 0.37 min) than at 20''C (13.6 ± 
0.38 min versus 17.3 ± 0.37 min). Evidence for effects of 
inbreeding on chill-coma recovery time were restricted 
to animals reared at 27°C (C: 25.5 ± 0.39 min = 12: 24.8 
± 0.39 min < II: 27.7 ± 0.39 min), while no effects of 
inbreeding were found in animals reared at 20°C (C: 14.4 
± 0.52 min = 12: 15.6 ± 0.55 min = II: 15.5 ± 0.52 min, 
Tukey HSD; significant inbreeding level * rearing tem- 
perature interaction). More specifically, the selection 
lines reared at 27°C caused the above pattern, while the 
unselected control lines showed neither at 20°C nor at 
27°C a significant response to inbreeding, and while the 
selection lines reared at 20°C showed an increase in 
chill-coma recovery time with increasing inbreeding 
level (significant selection regime * inbreeding level * 
rearing temperature interaction; Figure 2A). 

Likewise, inbreeding effects were only detectable in 
animals acclimated to 27''C (C: 24.8 ± 0.46 min = 12: 
24.9 ± 0.48 min < II: 27.6 ± 0.46 min), but not in those 
acclimated to 20°C (C: 15.1 + 0.46 min = 12: 15.5 ± 
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Figure 2 Variation in chill-coma recovery time (means -i- 1 SE) in Bicyclus anynana. Chill-coma recovery time in relation to rearing 
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0.48 min = II: 15.7 ± 0.46 min; Tukey HSD; significant 
inbreeding level * adult temperature interaction). The 
lack of a comparable response in the butterflies accli- 
mated to 20°C is caused by the unselected controls, 
showing an increase in recovery time with increasing 
inbreeding level (significant selection regime * inbreed- 
ing level * adult temperature interaction; Figure 2B). 
Overall, inbreeding effects were thus only detectable 
in the animals reared at and acclimated to 27°C, but not 
in any other rearing by adult temperature group (sig- 
nificant inbreeding level * rearing temperature * adult 
temperature interaction; Figure 2C). The significant 
adult temperature by sex interaction indicates that males 
tended to show longer recovery times than females at 
the adult temperature of 20°C (16.0 ± 0.38 min vs. 14.9 
± 0.38 min), but shorter ones than females at 27°C (25.2 
± 0.37 min vs. 26.3 ± 0.39 min). This pattern was 
restricted to animals reared at 27°C though, while sexual 
differences were absent in animals reared at 20°C (sig- 
nificant rearing temperature * adult temperature * sex 
interaction; Figure 2D). 



When analyzing the animals reared at 20°C separately, 
selection regime, adult temperature and replicate line 
significantly affected chill-coma recovery time, but not 
inbreeding level and sex (Table 3B). Butterflies from the 
selection lines recovered quicker compared to unse- 
lected control butterflies (13.2 ± 0.70 min < 17.2 ± 
0.66 min), and adults kept at 20°C recovered much faster 
than those kept at 27°C (10.1 ± 0.32 min < 20.3 ± 
0.33 min). All patterns are thus in full agreement with 
the above analysis. The interaction between selection re- 
gime and adult temperature indicates that differences in 
chill-coma recovery time between selection and control 
lines were more pronounced at the higher (17.4 ± 
0.47 min vs. 23.3 ± 0.45 min) compared to the lower 
adult temperature (9.0 ± 0.47 min vs. 11.1 ± 0.45 min). 

In the animals reared at 27°C, chill-coma recovery 
time was significantly affected by selection regime, adult 
temperature, and replicate line, but not by adult feeding 
regime, sex and inbreeding level (Table 3C). Butterflies 
from the selection lines had shorter recovery times than 
those from the unselected control lines (23.3 ± 0.99 min 
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Table 4 Results of nested ANOVAs for longevity 



A) Longevity 


DF 


MS 


F 


P 


Selection Regime 


1,6 


234.0 


0.6 


0.4616 


Replicate Line [Sel. & Inbr.] 


6,3744 


404.9 


2.8 


0.0109 


nbreeding Level 


2,6 


146.3 


04 


0.6967 


Rearing Temperature 


1 ,3744 


264262.0 


1 806.6 


<0.0001 


Sex 


1 ,3744 


20.7 


0.14 


0.7070 


Selection Reg. * Inbreeding Level 


2,6 


2020.7 


5.3 


0.0440 


Selection Reg. * Rearing Temp. 


1 ,3744 


93.4 


0.6 


04242 


Selection Reg. * Sex 


1 ,3744 


52.2 


04 


0.5505 


nbreeding Level * Rearing Temp. 


2,3744 


440.1 


3.0 


0.0495 


nbreeding Level * Sex 


2,3744 


153.7 


1.0 


0.3498 


Rearing Temp. * Sex 


1,3744 


5324.0 


36.4 


<0.0001 


jci. r\cy. iiiuiccu. rvccii. iciiiij. 


2 3744 


472 0 


3 2 


0.0398 


jci. r\cy. II luiccu. jca 


2 3744 


940 0 


54 


0.0016 


Sel. Reg. * Rearing Temp. * Sex 


1,3744 


1.5 


<0.1 


0.9202 


nbreed * Rear. Temp * Sex 


2,3744 


178.4 


1.2 


0.2954 


FrMir-\A/3\/ 1 ntpra("tif-in 
UUI Way IIILCId(^LIUII 


2 3744 


932 1 


54 


0.001 7 




3744 


146 3 






B) Animals reared at 20**C 


DF 


MS 


F 


P 


Selection Regime 


1,6 


287.5 


0.3 


0.6267 


nbreeding Level 


2,6 


174.0 


0.2 


0.8570 


Replicate Line [Sel. Reg. & Inbr.] 


6,1312 


1 1 06.3 


4.8 


<0.0001 


Adult Temperature 


1,1312 


34430.6 


148.6 


<0.0001 


Sex 


1,1312 


2568.9 


11.1 


0.0009 


Selection Reg. * Inbreeding Level 


2.6 


1712.1 


1.6 


0.2864 


Selection Reg. * Adult Temp. 


1,1312 


1 726.0 


7.5 


0.0064 


o/"tif-in Rfin * ^(^v 
JCICl^LIUI 1 r\trU . JCA 


11312 


7 3 


<0 1 


0 8589 


in hfpprl i nr~i 1 * AHi i t Ttimpi 
lllUICCUIIiy LCVCI AlUUIL ICIIILJ. 


2 1312 


3034 3 


1 3 5 




nbreeding Level * Sex 


2,1312 


131.9 


0.6 


0.5659 


Aril i t Tf=imnpr3ti irp * "^ov 
AAU U 1 L 1 CI 1 1 UCI a LU 1 C JCA 


11312 


2752 2 


1 2 6 


0.0004 


Rpn * lnnrf=ifiri * ^pv 
jci. ncy. II luiccru. JCA 


2 1312 


1 395 7 


50 


0.0025 




1312 


231 6 






C) Animals reared at 27°C 


DF 


MS 


F 


p 


Selection Regime 


1,6 


38.0 


0.1 


0.7603 


nbreeding Level 


2,6 


142.8 


04 


0.6974 


Replicate Line [Sel. Reg. & Inbr.] 


6,2394 


37. 


7.0 


<0.0001 


Adult Temperature 


1,2394 


1 177.7 


21.8 


<0.0001 


Adu t Feeding Regime 


1,2394 


52827.7 


9772 


<0.0001 


Sex 


1,2394 


3684.1 


68.1 


<0.0001 


Selection Reg. * Inbreeding Level 


2,6 


475.6 


1.3 


0.3451 


Selection Reg. * Adult Temp. 


1,2394 


104.6 


1.9 


0.1644 


Selection Reg. * Adult Feed. 


1,2394 


49.3 


0.9 


0.3399 


Selection Reg. * Sex 


1,2394 


15.4 


0.3 


0.5936 


Inbreeding Level * Adult Temp. 


2,2394 


37.1 


0.7 


0.5033 


nbreeding Level * Adult Feed. 


2,2394 


112.8 


2.1 


0.1244 



Table 4 Results of nested ANOVAs for longevity 

(Continued) 



Inbreedi 


ng Level * Sex 


2,2394 


137.9 


2.5 


0.0782 


Adult Temp. * Adult Feed. 


1,2394 


451.8 


84 


0.0039 


Adult Temp. * Sex 


1,2394 


114.8 


2.1 


0.1451 


Adult Feed. * Sex 


1,2394 


1481.2 


274 


<0.0001 


Sel. Reg 


* Inbreed. * Adult Temp. 


2,2394 


174 


0.3 


0.7248 


Sel. Reg 


* Inbreed. *Adult Feed. 


2,2394 


532.5 


9.9 


<0.0001 


Sel. Reg 


* Inbreed. * Sex 


22394 


0.7 


<0.1 


0.9879 


Sel. Reg 


* Adult Temp. * Sex 


2,2394 


1.7 


<0.1 


0.8610 


Sel. Reg 


* Adult Temp. * AdFeed. 


1,2394 


128.8 


24 


0.1229 


Sel. Reg 


* Adult Feed. * Sex 


1,2394 


0.7 


<0.1 


0.9110 


Inbreed. 


* Adult Temp. * AdFeed. 


2,2394 


144 


0.3 


0.7659 


Inbreed. 


* Adult Temp. * Sex 


2,2394 


161.9 


3.0 


0.0502 


Inbreed. 


* Adult Feed. * Sex 


2,2394 


164.2 


3.0 


0.0482 


Adult Temp. * Adult Feed. * Sex 


1,2394 


87.1 


1.5 


0.2044 


Inbreed. 


♦ AdT * AdF * Sex 


2,2394 


222.8 


4.1 


0.0163 


Error 




2394 


54.0 







Results for the effects of selection regime (selection vs. control), inbreeding 
level (C, II, 12}, replicate line, replicate cage, rearing temperature (20 vs. 27°C), 
and sex on longevity in Bicyclus anynana (A). In Tables B (animals reared at 
20°C) and C (animals reared at 27''C) data were analyzed separately for each 
rearing temperature, as different feeding treatments were employed in 
animals reared at 27°C only, and as different adult temperatures were 
employed in animals reared at 27°C only temporarily. Consequently, the 
analysis in A was restricted to the animals fed ad libitum and l<ept as adults at 
27°C, to mainly explore effects of rearing temperature. The factor rearing 
temperature was replaced by the factor adult temperature in analyses B and 
C, and the factor adult feeding regime (control vs. starvation) was added to 
Table C. For details see Material and Methods. Replicate line (random factor) 
was nested within selection regime and inbreeding level throughout. 
IVlinimum adequate models were constructed by sequentially removing 
non-significant interaction terms. Significant P-values are given in bold. 

VS. 28.8 ± 0.99 min), and individuals acclimated to the 
lower temperature recovered faster than those accli- 
mated to the higher temperature (20.9 ± 0.36 min < 31.2 
± 0.35 min). The significant interactions found here ba- 
sically reflect the patterns described above, and are thus 
not described here except for both interactions involving 
feeding regime. The response to adult feeding regime 
varied largely among inbreeding levels and sexes, with 
feeding regime having positive to negative effects (sig- 
nificant inbreeding level * adult feeding regime * sex 
interaction; Figure 2E). Furthermore, food stress had 
generally either little or positive effects on chill-coma re- 
covery times in the adult temperature by sex groups, 
while it slightly decreased cold stress resistance in the 
females acclimated to 27°C (significant adult temperature 
* adult feeding regime * sex interaction; Figure 2F). 

Longevity 

In the analysis addressing effects of rearing temperature, 
longevity was significantly affected by rearing 
temperature and replicate line, but not by selection re- 
gime, inbreeding level, and sex (Table 4A). Animals 
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reared at 20°C lived longer compared to animals reared 
at 27°C (29.9 ± 0.33 d > 12.3 ± 0.25 d). The significant 
selection regime by inbreeding level interaction indicates 
that longevity tended to be longer in the unselected con- 
trols than in the selection lines in the outbred controls 
(23.2 ± 0.78 d vs. 19.7 ± 0.81 d), but shorter in inbreed- 
ing group 1 (20.5 ± 0.81 d vs. 22.0 ± 0.79 d), while being 
very similar in inbreeding group 2 (20.6 ± 0.84 d vs. 
20.9 ± 0.87 d). These patterns, however, were largely 
restricted to animals reared at 20°C and to males (sig- 
nificant selection regime * inbreeding level * rearing 
temperature interaction and significant selection regime 
* inbreeding level * sex interaction; Figure 3A and 3B). 
Consequently, significant differences between selection 
and control lines were exclusively found in outbred 
males reared at 20°C (35.5 ± 1.1 d vs. 26.3 ± 1.1 d; sig- 
nificant four- way interaction). When being reared at 
20°C, longevity was longest in inbreeding 1 group 
(C: 30.0 ± 0.56; II: 30.7 ±0.56 d; 12: 29.0 ± 0.61 d), while 
it was shortest in this group when being reared at 27°C 
(C: 12.8 ± 0.42; II: 11.7 ± 0.42 d; 12; 12.4 ± 0.43 d; sig- 
nificant inbreeding level * rearing temperature inter- 
action). Sex differences were only found in animals reared 
at 20°C (males: 31.1 ± 0.47 d > females: 28.8 ± 0.47 d), 
while there was no significant sex difference in animals 
reared at 27°C (males: 11.0 ± 0.35 d = females: 13.7 ± 
0.35 d; Tukey HSD; significant rearing temperature * sex 
interaction). 

When considering the animals reared at 20°C separ- 
ately, adult temperature, sex, and replicate line signifi- 
cantly affected longevity, but not selection regime and 



inbreeding level (Table 4B). Adults acclimated to 20°C 
lived longer than adults acclimated to 27°C (34.7 ± 0.57 
d > 24.4 ± 0.62 d), and males lived longer than females 
(30.9 ± 0.60 d > 28.1 ± 0.60 d). The sex difference was 
restricted to the adult temperature of 27°C though 
(males: 27.3 ± 0.87 d > females: 21.5 ± 0.88 d), while 
there was no significant difference at 20°C (males: 34.6 ± 
0.82 d = females: 34.7 ± 0.80 d; Tukey HSD; signifi- 
cant adult temperature * sex interaction). The signifi- 
cant interaction between selection regime and adult 
temperature indicates that longevity was longer in unse- 
lected controls compared to selection lines at the lower 
adult temperature (36.2 ± 0.80 d vs. 33.1 ± 0.82 d), but 
vice versa at the higher adult temperature (23.8 + 0.86 d 
vs. 24.9 ± 0.88 d). Longevity decreased with increasing 
inbreeding level when animals were acclimated to 20°C 
(C: 37.2 ± 0.98 d = II: 35.8 ± 0.94 d > 12: 31.0 ± 1.05 d), 
but tended to increase with increasing inbreeding level 
when animals were acclimated to 27°C (C: 22.2 ± 1.0 d = 
II: 24.4 ± 1.06 d = 12: 26.6 ± 1.13 d; Tukey HSD; signifi- 
cant inbreeding level * adult temperature interaction). 
Finally, males of the outbred control and the inbreeding 
1 group showed quite substantial though opposing dif- 
ferences between selection lines and unselected controls, 
while differences where much smaller in other groups 
(significant selection regime * inbreeding level * sex 
interaction; Figure 3C). 

In the animals reared at 27°C, longevity was signifi- 
cantly affected by adult temperature, adult feeding re- 
gime, sex and replicate line, but not by selection regime 
and inbreeding level (Table 4C). Adults acclimated to 
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20°C lived longer than adults acclimated to 27°C (13.5 ± 
0.21 d > 12.1 ± 0.21 d), control animals lived longer 
compared to food-stressed animals (17.5 ± 0.22 d > 8.1 
± 0.20 d), and females lived longer than males (14.1 ± 
0.21 d > 11.6 ± 0.21 d). Effects of adult temperature were 
more pronounced in control (20°C: 18.7 ± 0.32 d vs. 27°C: 
16.4 ± 0.31 d) than in food-stressed animals (20°C: 
8.4 ± 0.29 d vs. 27°C: 7.8 ± 0.28 d; significant adult 
temperature * adult feeding regime interaction). Sex dif- 
ferences, in contrast, were larger in food-stressed 
(female: 10.2 ± 0.29 d vs. male: 6.1 ± 0.28 d) than in 
control animals (female: 18.0 ± 0.32 d vs. male: 17.1 ± 
0.31 d; significant adult feeding regime * sex inter- 
action). While in food-stressed animals longevity was 
generally very similar across selection regimes and 
inbreeding levels, the feeding control groups showed di- 
vergent patterns between unselected controls and selec- 
tion lines (Figure 3D). While in the latter groups 
longevity increased with increasing inbreeding level in the 
selection lines (C: 16.1 ± 0.56 d = U: 16.9 ± 0.56 d < 12: 
18.8 ± 0.54 d), it decreased with increasing inbreeding 
level in the unselected control lines (C: 20.2 ± 0.51 d > II: 
16.4 ± 0.56 d = 12: 16.8 ± 0.56 d; significant selection 
regime * inbreeding level * adult feeding regime inter- 
action). The significant inbreeding level * adult feeding 
regime * sex interaction shows that inbreeding effects 
were only present in food control males, in which indi- 
viduals from inbreeding level 1 lived shortest (II: 15.5 ± 
0.56 d < 12: 17.0 ± 0.55 d = C: 18.7 ± 0.53 d). Variation 
between inbreeding groups though was restricted to food 
control males kept at 20°C (II: 16.1 ± 0.79 d; 12: 18.1 ± 
0.78 d; C: 21.6 ± 0.76 d; significant inbreeding level * adult 
temperature * adult feeding regime * sex interaction). 

Discussion 

Variation In developmental traits 

As predicted (Table 1), selection regime neither affected 
larval time, pupal time, larval growth rate or pupal mass, 
thus suggesting that genetic variation in adult cold toler- 
ance operates largely independent of these developmen- 
tal traits in B. anynana (see [17]; cf [16,42]). However, 
selection regime was involved in a total of six interac- 
tions with other factors (larval time: 3; pupal time: 1; 
growth rate: 1; pupal mass: 1; Table 2). The patterns 
revealed through these interactions though do not fun- 
damentally change the conclusion drawn above, with 
perhaps the following exception. Larval development 
time was slightly longer in control than in selection lines 
at the lower rearing temperature only. This may indicate 
that the lines selected for increased adult cold toler- 
ance may be slightly better adapted to deal with devel- 
opmental temperatures below the optimal range. Note 
though that this pattern was not consistent across all 
treatment groups (significant 3-way interaction). All 



other interactions indicated inconsistent, group-specific 
responses and did therefore not reveal any straightfor- 
ward effects of selection regime on developmental 
traits. 

Likewise, inbreeding yielded only marginal effects, sig- 
nificantly affecting pupal mass only. Pupal mass was 
highest in inbreeding group 1, intermediate in the out- 
bred control, and lowest in inbreeding group 2. Directly 
after the full-sib matings, in contrast, larval development 
time, larval growth rate, and pupal mass showed 
inbreeding depression [17]. The fact that much of these 
initial differences disappeared indicates fitness rebounds 
caused by subsequent random mating (for 14 genera- 
tions: 10 during the selection experiments and 4 after- 
wards) as was expected [17]. Additive theory would 
predict a complete rebound under such circumstances 
[43,44], such that the pattern found for pupal mass may 
reflect a chance effect caused by random genetic drift ra- 
ther than inbreeding depression. Note in this context 
that the 11 interactions in which inbreeding was 
involved in also indicated erratic variation rather than 
any systematic patterns. 

In contrast to the above results, sexual differences and 
variation induced by different rearing temperatures resem- 
ble predicted and well-known patterns for B. anynana and 
beyond [45,46]. Specifically, lower temperatures resulted 
in increased development time and pupal mass, but in 
decreased larval growth rates. This is consistent with the 
temperature-size rule and the in general temperature- 
dependent growth of ectotherms [46-48] . Males compared 
to females showed shorter larval but longer pupal times 
(though with the latter not affecting earlier male emer- 
gence; protandry selection), higher larval growth rate, and 
lower pupal mass (fecundity selection in females). Sex dif- 
ferences were more pronounced at 20°C (larval time) or at 
27°C (pupal time, larval growth rates). While the former 
has been reported previously and likely reflects larger ab- 
solute differences owing to the overall longer development 
time at lower temperatures [49], we have no explanation 
for the latter. 

Variation in chill-coma recovery time 

Cold tolerance was strongly affected by selection regime, 
with recovery times being overall by 21% shorter in the 
selection than in the control lines. Directly after the 
course of selection the difference between selection 
regimes was 29% [17], indicating an ongoing conver- 
gence towards control line levels. This suggests that the 
selected genes and alleles had not yet become fixed in 
the selection lines, and furthermore a lack of selection 
for maintaining increased cold tolerance under labora- 
tory conditions. Interestingly, the differences between 
selection and control lines tended to be more pro- 
nounced when animals had been reared at 20°C (23%) as 
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compared to 27°C (18%; Table 3A). This suggests that 
the selection lines might show a somewhat stronger re- 
sponse to rearing temperature (by 43% shorter chill- 
coma recovery time at 20°C compared to 27°C) than the 
control lines (by 40% shorter). 

Regarding interactive effects, selection regime was 
involved in a total of five significant interactions. While 
three of these did not reveal any conclusive patterns 
regarding the role of selection regime, the remaining 
two interactions suggest that the differences between se- 
lection and control lines were larger at the higher as 
compared to the lower adult temperature (Table 3 A and 
B). However, while this is true in absolute terms 
(5.8 minutes at 27°C vs. 3.7 minutes at 20°C for all indi- 
viduals; 5.9 minutes at 27°C vs. 2.1 minutes at 20°C for 
the individuals reared at 20°C), larger relative differences 
were found in the animals reared at 20°C only (27°C: 
20% and 20°C: 21% for all; 27°C: 25% and 20°C: 19% for 
individuals reared at 20°C). Concomitantly plastic 
responses are quite similar across selection regimes, with 
the selection lines showing by 41% (all) and 48% (reared 
at 20°C) shorter chill-coma recovery times at 20°C com- 
pared to 27°C, and the control lines by 40% and 52% 
shorter ones. We therefore conclude that the extent of 
the plastic responses in cold tolerance is comparably 
high across both selection regimes, thus challenging the 
main hypothesis we wanted to test here. 

Effects of different rearing and adult temperatures 
were strong and persistent, with animals reared or accli- 
mated at the lower temperature showing increased cold 
tolerance compared to those kept at the higher 
temperature (cf Table 1). Similar temperature-induced 
changes in cold tolerance have been documented in a 
large variety of insects, suggesting a universal pattern of 
adaptive phenotypic plasticity [23,50]. In contrast, no 
consistent effects of adult feeding treatment on cold tol- 
erance were found, despite large variation among 
inbreeding levels, sexes, and temperature groups. This is 
in line with our predictions based on earlier results [23], 
which indicates that cold tolerance is at least not 
strongly limited by adult food intake. Furthermore, our 
results revealed no sex differences in cold tolerance as 
was expected. Likewise, testing for sex differences in 
temperature stress resistance in an earlier study yielded 
a non-significant result in 18 out of 22 statistical ana- 
lyses, suggesting that, in general, both sexes are equally 
tolerant in B. anynana [23]. Except from some group- 
specific responses, no effects of (earlier) inbreeding on 
cold tolerance were found. 

Variation In longevity 

Longevity was not affected by selection regime, indicat- 
ing that there is no trade-off between increased cold 
stress resistance and life span in B. anynana (cf [34]). 



This notion is further supported by the erratic patterns 
indicated by the seven interactions in which selection re- 
gime was involved in. Studies on Drosophila also yielded 
inconclusive results regarding genetic correlations be- 
tween cold tolerance and longevity (e.g. [2,15,16,42]). 
Temperature effects, in contrast, were straight-forward 
and as predicted. Animals reared at or acclimated to 20°C 
lived generally longer than animals reared or kept at the 
higher temperature, thus indicating accelerated rates of 
ageing at warmer temperatures owing to higher metabolic 
rates [2,47,51]. Adult food stress had as expected a large 
impact on longevity, with a lack of food substantially re- 
ducing life span [36,37]. Not surprisingly, differences in 
longevity between acclimation temperature groups were 
more pronounced in control compared to food-stressed 
animals, as all latter individuals died within relatively 
short time spans. Once again, inbreeding yielded no con- 
clusive effects. 

Regarding sex differences, females lived longer than 
males when being reared at 27°C, as has been found in 
several earlier studies [17,38]. The sex difference was 
more pronounced in food-stressed than in control ani- 
mals, which may suggest a re-allocation of resources 
from reproduction to somatic maintenance under food 
stress in females. When being reared at 20°C, however, 
males lived longer than females, but only at the higher 
adult temperature. This finding may suggest that the re- 
spective temperature regime does not only affect longev- 
ity as such, but influences sexes in different ways. 

Conclusions 

Regarding effects of temperature, feeding treatment, sex, 
inbreeding and selection regime our results are in good 
agreement with the hypotheses outlined in Table 1 and 
thus with earlier findings obtained in B. anynana and 
other insects. However, our principal intention here was 
to explore plastic responses in populations of different 
genetic background and thus genotype-environment 
interactions. Regarding this matter we have indeed 
detected a number of interactions in which selection re- 
gime was involved in, but most of these did not reveal 
consistent patterns. We have found some evidence 
though that the lines selected for increased cold toler- 
ance as compared to control lines (1) developed slightly 
faster at suboptimally low temperatures, (2) showed lar- 
ger differences in cold tolerance compared to controls 
when being reared at the lower temperature, (3) and 
showed more pronounced differences in cold tolerance 
compared to controls at an adult temperature of 27°C. 
The former findings may indicate that the lines selected 
for increased adult cold tolerance may be slightly better 
equipped to deal with lower temperatures, such that any 
advantages become exaggerated when experiencing sub- 
optimally low temperatures during development. While 
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thus the plastic response to different rearing tempera- 
tures seemed to be somewhat larger in the selection 
compared to the control lines, plastic responses to dif- 
ferent adult temperatures were very similar across selec- 
tion regimes. Our results thus suggest that plastic 
capacities were, if anything, only marginally affected by 
directional selection, and that plasticity was largely inde- 
pendent of the respective trait mean. This suggests that 
the genetic architecture of the trait mean operates 
largely independently of the architecture underlying 
plastic responses. Thus, we found no evidence for an 
interference of genetic adaptation with plastic capacities 
(but see [52] for an interspecific comparison). 

However, one should keep in mind that our study 
draws on variation across selection hnes, which had been 
established from a single laboratory source population. 
While using a single gene pool to experimentally estab- 
lish genetically differentiated populations has clear 
advantages, it also involves some limitations which 
should be acknowledged. Notably, we do not know 
whether our results are transferable to other populations 
(and species), and whether they might be affected by the 
specific population history (e.g. inbreeding, founder 
effects, laboratory adaptation). Several traits have been 
shown to readily respond to laboratory conditions, 
though essentially nothing is known about the effects of 
laboratory adaptation on plastic responses [53,54]. For 
logistic reasons such limitations apply to many labora- 
tory experiments. At least some studies indicate that la- 
boratory results may nevertheless have clear relevance 
for field conditions [54]. In any case our results clearly 
show that there is not necessarily a link between trait 
value and the magnitude of plastic responses. 
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